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2} ■ Abstract 

■ The fragmentation function of vector particle into possible bound S-wave states with a 

heavy antiquark is calculated in the leading order of perturbative QCD for the high energy 
^ I processes at large transverse momenta with the different behaviour of anomalous chromo- 

<^ ' magnetic moment. One- loop equations are derived for the evolution of fragmentation function 

moments, which is caused by the emission of hard gluons by the vector particle. The integrated 
probabilities of fragmentation are given. The distribution of bound state over the transverse 
momentum with respect to the axis of fragmentation is calculated in the scaling limit. 

(N ■ 1 Introduction 
cn 

00 

I An interesting problem concerning properties of interaction beyond the Standard Model is to study 

^ ' production of hadrons containing leptoquarks [|l| which are scalar and vector color-triplet particles 

appearing in Grand Unification Theories, provided that their total width is much less than the 
QCD-confinement scale, Tlq -C Aqcd- Recently the production of (gLQ)-baryons in the case of 
p . scalar leptoquark was discussed 0. 

I In this work we study the high energy production of heavy leptoquarkonium containing a vector 

' leptoquark. These results can be used to calculate the fragmentation of vector local diquarks into 

baryons (a similar approach was applied to the production of Qccc in 0)- For the sake of convenience 
^ . the local color-triplet vector field will be referred to as the leptoquark in this paper. 

New problem arising in this case is a choice of the lagrangian for the vector leptoquark interaction 
with gluons. Indeed, to the lagrangian of a free vector field 

where Hf^^, = d^Uy — dpU^^, t/^ is the vector complex field with derivatives substituted by covariant 
ones, we can add the gauge invariant term proportional to 

where S"^/^ = 1/2 (5° 5^ — 5" 5^) is the tensor of spin, G^" is the gluon field strength tensor. It 
leads to the appearance of a parameter in the gluon-leptoquark vertex (the so-called anomalous 
chromomagnetic moment, see Section 2). In this work we discuss the production of a 1/2-spin 
bound state containing the heavy vector particle at various values of this parameter. 
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At high transverse momenta, the dominant production mechanism for the heavy leptoquarko- 
nium bound states is the leptoquark fragmentation, which can be calculated in perturbative QCD 
[^] after the isolation of soft binding factor extracted from the non-relativistic potential models 
[§, ^. The corresponding fragmentation function is universal for any high energy process for the 
direct production of leptoquarkonia. 

In the leading a^-order, the fragmentation function has a scaling form, which is the initial con- 
dition for the perturbative QCD evolution caused by the emission of hard gluons by the leptoquark 
before the hadronization. The corresponding splitting function differs from that for the heavy quark 
because of the spin structure of gluon coupling to the leptoquark. 

In this work the scaling fragmentation function is calculated in Section 2 for two different cases 
of the anomalous chromomagnetic moment behaviour. The limit of infinitely heavy leptoquark, 
'TTiLQ ^ oo, is obtained from the full QCD consideration for the fragmentation. The distribution of 
bound state over the transverse momentum with respect to the axis of fragmentation is calculated in 
the scaling limit in Section 3. The splitting kernel of the DGALAP-evolution is derived in Section 4, 
where the one-loop equations of renormalization group for the moments of fragmentation function 
are obtained and solved. These equations are universal, since they do not depend on whether 
the leptoquark will bound or free at low virtualities, where the perturbative evolution stops. The 
integrated probabilities of leptoquark fragmentation into the heavy leptoquarkonia are evaluated in 
Section 5. The results are summarized in Conclusion. 



2 Fragmentation function in leading order 

The contribution of fragmentation into direct production of heavy leptoquarkonium has the form 

da[lH{p)]= / dz dd-[LQ{p/z),jj\ Dlq^Ih{z,IJ^), 
Jo 

where da is the differential cross-section for the production of leptoquarkonium with the 4-momen- 
tum p, da is that of the hard production of leptoquark with the scaled momentum p/z, and D is 
interpreted as the fragmentation function depending on the fraction of momentum carried out by 
the bound state. The value of /i determines the factorization scale. In accordance with the general 
DGLAP-evolution, the /i-dependent fragmentation function satisfies the equation 



dDLQ^l„{z,fl) 



J z y 



where P is the kernel caused by emission of hard gluons by the leptoquark before the production of 
heavy quark pair. Therefore, the initial form of fragmentation function is determined by the diagram 
shown in Fig.|l], and, hence, the corresponding initial factorization scale is equal to /i = 2mQ. 
Furthermore, this function can be calculated as an expansion in as{2mQ). The leading order 
contribution is evaluated in this Section. 

Consider the fragmentation diagram in the system, where the momentum of initial leptoquark 
has the form q = (go, 0, 0, qs) and the leptoquarkonium one is p, so that 

g2 = pi = M^. 
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Figure 1: The diagram of leptoquark fragmentation into the heavy leptoquarkonium. 

In the static approximation for the bound state of leptoquark and heavy quark, the quark mass is 
expressed as uiq = rM, and the leptoquark mass equals m — {1 — r)M — fM. The gluon-vector 
leptoquark vertex has the form 

= -igs^lQuuiq + fp)a - Qf^aiii + se)rp - eeq)^ - + se)q - eerp)^], (2) 

where ae is the anamalous chromomagnetic moment, t"' is the QCD group generator. The sum over 
the vector leptoquark polarizations with the momentum q (g^ = s) depends on the choice of the 
gauge of free field lagrangian (for example, the Stueckelberg gauge) , but the fragmentation function 
is a physical quantity and has not to depend on the gauge parameter changing the contribution of 
longitudinal components of the vector field. So, the sum over polarization can be taken in the form 

s 

The matrix element of the fragmentation into the baryon with the spin of 1/2 has the form 



M = - ''l^: rf{f-m-f n<lUP{fp),,TZ^p^, - Mh^^H, Ml (3) 
where the sum over the gluon polarization is written down in the axial gauge with n — (1, 0, 0, —1) 



k ■ n 



with k — g — (1 — r)p. The spinors Ih and q correspond to the leptoquarkonium and heavy 
quark associated to the fragmentation. A4q denotes the matrix element for the hard production 
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of leptoquark at high energy, -R(O) is the radial wave-function at the origin. The matrix element 
squared and summed over the helicities of particles in the final state has the form 



,-,,2 



|M| =W^,M^M^ 

In the limit of high energies q ■ n —>■ oo, W^^ behaves like 

W^, = -g^,W + R^,, (4) 

where R^y can depend on the gauge parameters and leads to scalar formfactor terms which are 

small in comparison with W in the limit of g ■ n — oo. Define 

p ■ n 

z = . 

q ■ n 

The fragmentation function is determined by the expression [0 

where W is defined in (^). The integral in the expression for the fragmentation function diverges 
logarithmically at a constant value of anamalous chromomagnetic moment if « does not equal — 1. 
In this work we consider two sets for the behaviour of anamalous chromomagnetic moment. 

Set I is defined by se = —1. Here we observe that the obtained fragmentation function coincides 
with that for the scalar leptoquark |^ except the factor of 1/3 

= SS^^ ■ {3 + 3.^ -(6 -10. + 2.= + 2.').+ 

2437r M-^r'^r^ [1 — rz^ I 

+ (3-10r + 14r2-10r^ + 3r^);z^)|, (5) 

which tends to 

24377 y° rriQ r I J 

at r ^ and y = (1 — (1— r)^) / (rz). The limit of D{y) is in agreement with the general consideration 
of 1/m-expansion for the fragmentation function where 

D{y) = -a{y) + b{y). 
r 

The dependence on y in a{y) is the same as for the fragmentation of the heavy quark into the 
quarkonium |7]]. 

Set II: se behaves like —1 + AM'^/{s — ^\q)- The obtained fragmentation function is 

Sag |/?(0)P z''{l-zf 
2437r 16M3r2f2 [1 - fzf 

|l6(3 + 3r2 -Qz+ lOrz - 2r'^z - 2r^z + 

+3^2 - 10;zV + 14A^ - lO^V^ + ?,z^r*) + 
+A(3A + 24r - QzA - IrzA - 'ilrz - IQr^z + 3z^A + 

+2z^rA + 8zV + 3zVM - 32z^r^ + 24zV^) | . (7) 



D(z) 



^In paper |2| an arithmetic error was made, which shghtly affects upon the finai result at small 
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D/N 




z 



Figure 2: The fragmentation function of leptoquark into the heavy leptoquarkonium, the A^-factor 
is determined by = ^gi^ljzTp-, the fragmentation function for Set I is shown by the dashed 
hne, the fragmentation function for Set II [A = 3) is given by the solid line at r = 0.02. 

It tends to 

8al |i?(0)|Ml/-l)' 



Diy) 



2437r 16y^ m| r 

1 16(8 + % + 3y^) + A{8A - 8yA + 3y^A - 64 + 16y)|, (8) 



at r — s> and y = {1 — {1 — r)z) / (rz). The perturbative functions in the leading Os-order are shown 
in Fig.0 at r = 0.02. We see hard distributions, which become softer with the evolution (see ref.0). 

3 Transverse momentum of leptoquarkonium 

In the system with an infinite momentum of the fragmentating leptoquark its invariant mass is 
expressed by the fraction of longitudinal leptoquark momentum z and transverse momentum with 
respect to the fragmentation axis pr (see Fig. |l]) as 

^ " + zil-z) " " + 
where t = px/M. The calculation of diagram in Fig. |T] gives the double distribution 

ds dz ' ' 
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where the Set I function T> has the form 

256a2 |i?(0)|2 m3 



For Set II, we find 

V{z,s) 



2 2 o 



|i?(0)|2 M3 



Slvr r^r^ [1 — rz]^(s — m^)'^ 
|8r(A - 4 + 4r)^(l - z + rz)^ + 
+2(-y4 - 4 - 4r + + 4z - rzA + 16r2; - 4r^z) 

(1 - z + rz)(A - 4 + 4r)^^ - 32(1 - ^)^(^^)'}- (10) 

The distribution of leptoquarkonium over the transverse momentum can be obtained by the inte- 
gration over z 

D(t) = / dzV{z,s) 







zil-zY 



For Set I we have 

m 



64a2 1^(0) p 



8l7r 3(1 -r)5M3 

|t(30r^ - 30r^ - 61tV + 45r2t2 + ssrH'^ _ 

(30r^ - 99rH^ - 54r^t^ + 27r^t2 + 9t^ + 18rt^ - 6rY + 



2A{2rH + rt' + rH') In ^'^^J^^f ^ ) }. (11) 

The distribution for Set II is given in Appendix. The typical form of distribution over the transverse 
momentum is shown in Fig. ^ 

4 Hard gluon emission 

The one-loop contribution can be calculated in the way described in the previous Sections. This 
term does not depend on the part of leptoquark-gluon vertex with the anomalous chromomagnetic 
moment, therefore the splitting kernel coincides with that for the scalar leptoquark. It equals 
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Figure 3: The distributions over the transverse momentum with respect to the axis of fragmentation, 
A^r factor is determined by Nt = |^ jj^^-;:y- The dashed hne represents the result for Set I, the 
sohd hne does it for Set II. 



where the "plus" denotes the ordinary action: dxf+{x) ■ g{x) = dxf{x) ■ [g{x) 
scalar leptoquark splitting function can be compared with that of the heavy quark 



gil)]. The 



Q- 



[x,n) 



l + x^ 



L 1 



X 



which has the same normalization factor at a; — > 1. Furthermore, multiplying the evolution equation 
by and integrating over z, one can get from (|^) the /i-dependence of moments a^n) up to the 
one-loop accuracy of renormalization group 



da 



(n) 



dinfj, 



Svr 



+ 



Hn), 



n> 1. 



(13) 



L2 n+1. 

At n = the right hand side of (0) equals zero, which means that the integrated probability of 
leptoquark fragmentation into the heavy leptoquarkonium is not changed during the evolution, and 
it is determined by the initial fragmentation function calculated perturbatively in Section 2. 
The solution of equation (p!3D has the form 



Hn) 



a 



3/3o 



n+1 



(14) 



-a^(yUo)J 

where one has used the one-loop expression for the QCD coupling constant 

"^^^^ (3oH^^/AQCDy 

where (3o = 11— 2n//3, with being the number of quark flavors with niq < /i < rriLQ- Relation (|1^) 
is universal one, since it is independent of whether the leptoquark is free or bound at the virtualities 
less than yUg- We can use the evolution for the fragmentation into the heavy leptoquarkonium. The 
leptoquark can lose about 20% of its momentum before the hadronization p|. 
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5 Integrated probabilities of fragmentation 



As has been mentioned above, the evolution conserves the integrated probabihty of fragmentation 
which can be calculated explicitly 

/•^ ^, , |i?(0)|2 

/ dz D(z) = — - J w(r). (15) 
J ^ ' Slvr 16m| ^ ^ ^ ^ 

For Set I, we have 

16[(8 + 15r - 60r2 + lOOr^ - 60r^ - 3r^) + 30r(l - r + + r^) In rl , , 

For Set II, we find {A = 3) 

w(r) = I (143 + 701r - 1882r^ + 3250r^ - 3245r^ + 2017r^ - 936r'^ - 48r'^) + 

+30r(25 - 21r + 43r^ + + 8r^ + 16r^) Inr)^^ — -. (17) 

J 15(1 — r)^ 

The w(r) functions are shown in Fig. ^ at low r. 

6 Conclusion 

In this work the dominant mechanism for the production of bound states of spin 1/2 of a local 
color-triplet vector field with a heavy antiquark is considered for high energy processes at large 
transverse momenta, where the fragmentation contributes as the leading term. We investigate 
two sets of the anamalous chromomagnetic moment behaviour. Set I is defined hj se = —1 (the 
expression for the fragmentation function diverges logarithmically at a constant value of anamalous 
chromomagnetic moment if se is not equal to —1). Here we observe that the obtained fragmentation 
function coincides with that for the scalar leptoquark up to the factor of 1/3 

^^'^ ~ 2437r M3rV-2 [1 - fzf V + '^^ ^^"^ ^ +Zr )z + 



-(3 - lOr + 14r^ - lOr^ + 3r^)z^) | , (18) 



D z) 



where r is the ratio of heavy quark mass to the mass of the bound state. In the infinitely heavy 
leptoquark limit, D{z) has the form, which agrees with what expected from general consideration 
of 1/m-expansion for the fragmentation function. Set II is defined by ae = — 1 + AM'^/[s — m\Q). 
The fragmentation function for Set II differs from that for Set I: 

8ag |i?(0)P z^{l-zY 
2437r 16M3r2f2 [1 - fzf 

|l6(3 + 3r2 -Qz+ lOrz - 2r'^z - 2r^z + 

+3^^ - lO^V + 14zV2 - lO^V^ + ?,z^r*) + 
+A{?,A + 24r - QzA - 2rzA - ?,2rz - IQr^z + 3z'^A + 

+2z^rA + 8zV + 3zVM - 32z^r^ + 24zV^)|. (19) 
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w(r) 




Figure 4: The w-functions for the leptoquark fragmentation into the heavy leptoquarkonium versus 
the fraction r = ttiq/M. The curves correspond to Sets as in Fig. |^. 

The distribution of bound state over the transverse momentum with respect to the axis of fragmen- 
tation is calculated in the scaling limit. The corresponding distribution functions are given by the 
expression (11) for Set I and the expression from Appendix. The hard gluon corrections caused by 
the splitting of vector leptoquark are taken into account so that the evolution kernel has the form 



4as(/i) 



2x 



X 



x,/i) = ^^ , (20) 



which results in the corresponding one-loop equation for the moments of fragmentation functions 
(see eqs.(13), (14)). The numerical estimates show that the probabilities of fragmentation into 
bound states with c- and b-quarks of heavy vector leptoquark with the mass about 400 GeV are 
of the order of lO"^^"^-*. This suppression makes the experimental observation of such states rather 
difficult. However we can use the perturbative expressions for the model of fragmentations into 
doubly heavy baryons, where the integrated probabilities are of the order of 10^*^^^^-*. The results 
will be used to investigate the fragmentation of doubly heavy vector diquarks into baryons elsewhere. 
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7 Appendix 

The expression for the distribution over the transverse momentum for Set II (se = —1 + '^^l ) is 
given by 

^ ' 8l7r 3(1 -r)7M3 

jsOr^t - 270r^t + 720r^t - 480r^t - 727tV + 2AAlrh^ - 

-2041r=^t3 - 65r^t^ + mArH^ - 272rH^ + 129t^ - 285tV + 
+4:77rH^ - mrH^ + 312rH^ - lUrH^ - (30r^ - 240r5 + 480r^ + 
+t^{-909r^ + 2412r^ - 45r^ - 1080r^ + 432r^) + 
+t^(171 + 288r - 492r2 + 696r=^ - 165r^ + 264r^ + 48r^) + 

+t^(129 - 156r + 321r2 - 168r3 + 144r^))arctgf ~ ^j^ ) + 

V r + / 

-M92tr(l - r)(-2t2(l - r^) + (1 - 4r)r2)ln(^^^^^:^^) |. (21) 
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